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Abstract 



Let Hi, H.2 be finite dimensional complex Hilbert spaces describing the states of 
two finite level quantum systems. Suppose pi is a state in Hi, i = 1,2. Let C(pi, p-i) 



be the convex set of all states p in Ti = H\®H2 whose marginal states in TL\ and H2 
are pi and p2 respectively. Here we present a necessary and sufficient criterion for a 
p in C{p\, P2) to be an extreme point. Such a condition implies, in particular, that 
for a state p to be an extreme point of C{pi , P2) it is necessary that the rank of p does 
not exceed (d\ + d\ — l) 5 , where di = dim Hi, i = 1,2. When Hi and H2 coincide 
with the 1-qubit Hilbert space C 2 with its standard orthonormal basis {|0 >, |1 >} 
and pi = P2 = \l it turns out that a state p £ C(^I, ^1) is extremal if and only if p 
is of the form |0 >< Q\ where \tt >= (|0 > \ip > +|1 > |V>i >) , {l^o >, >} 
being an arbitrary orthonormal basis of C 2 . In particular, the extremal states are 
the maximally entangled states. 
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1 Introduction 



One of the well-known problems of classical probability theory is the determination 
of the set of all extreme points in the convex set of all probability distributions in 
a product Borel space (X x Y, T x Q) with fixed marginal distributions p and v on 
(X, T) and (Y : Q) respectively. Denote this convex set by C(p,v). When X = Y = 
{1, 2, . . . , n} , T — Q is the field of all subsets of X and p — v is the uniform distribution 
then the problem is answered by the famous theorem of Birkhoff [1] that the set of 
extreme points of the convex set of all doubly stochastic matrices of order n is the set 
of all permutation matrices of order n. Problems of this kind have a natural analogue in 
quantum probability. Suppose Hi and H 2 are finite dimensional complex Hilbert spaces 
describing the states of two finite level quantum systems Si and S 2 respectively. Then 
the Hilbert space of the coupled system is Hi <8> H 2 . Suppose pi is a state of Si in 
Hi, i = 1,2. Any state p in Si 2 yields marginal states Tr Wa p in Hi and Tr Wl p in H2 where 
Tr w . is the relative trace over Hi. Denote by C (pi, p 2 ) the convex set of all states p of the 
coupled system Si 2 whose marginal states in Hi and H 2 are pi and p 2 respectively. One 
would like to have a complete description of the set of all extreme points of C (pi,p 2 ) . 
In this paper we shall present a necessary and sufficient criterion for an element p in 
C (pi, p 2 ) to be an extreme point. This leads to an interesting (and perhaps surprising) 
upper bound on the rank of such an extremal state p. Indeed, if p is an extreme point of 
C (pi, p 2 ) then the rank of p cannot exceed (df + d\ — l) 2 where di = dim Hi. Note that 
the rank of an arbitrary state in Hi ®H 2 can vary from 1 to did 2 . When Hi = H 2 = C 2 , 
{|0 >, |1 >} is the standard (computational) basis of C 2 and pi — p 2 — |J it turns out 
that a state p in C (|/, is extremal if and only if p has the form \Q >< Q\ where 
\Q >= -j= (|0 > l^o > +|1 > IV'i >) ! {I^o >j IV'i >} being any orthonormal basis of C 2 . 
These are the well-known maximally entangled states. 

This work was done by the author during his visit to the University of Greifswald 
during 17 June - 16 July under a DST (India) - DAAD (Germany) project between the 
Indian Statistical Institute and the mathematics department of the University of Greif- 
swald. The author is grateful to these organisations for their generous support. The 
hospitality extended by the colleagues of the Quantum Probability group in the Univer- 
sity of Greifswald and, particularly, Michael Schurmann is gratefully acknowledged. 
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2 Extreme points of the convex set C (pi, P2) 



In the analysis of extreme points in a compact convex set of positive definite matrices 
the following proposition plays an important role [5]. See also [2-4]. 



Proposition 2.1 Let p be any positive definite matrix of order n and rank k < n. Then 
there exists a permutation matrix a of order n, a k x (n — k) matrix A and a strictly 
positive definite matrix K of order k such that 



apa 



1 _ 


K 


KA 




A^K 


A^KA 



(2.1) 



If, in addition, p = \{p' + p") where p' and p" are also positive definite matrices then 
there exist positive definite matrices K', K" of order k such that 



4t - 

ap^a 



1 


K* 


K*A 




A^K* 


A^K*A 



(2.2) 



where # indicates / and //. 



Proof: Choose vectors Uj G C n , i — 1, 2, . . . , n such that 

p= (((uiluj))), i,j G {1,2, ...,n}. 

Since rank p = k, the linear span of all the u^s has dimension k. Hence modulo a per- 
mutation a of {1, 2, . . . , n} we may assume that U\, 112, ■ ■ ■ , Uk are linearly independent 
and 

u k +j — aijUi + a 2 jU 2 H V a kj u k , l<j<n — k. (2.3) 

Putting 

K = ({{ Ui \ Uj ))), i.j E 1,2,..., k, 

A = ((a i:j )), i = 1,2,..., k; j = l,2,...,n-k 

and denoting by the same letter a, the permutation unitary matrix of order n corre- 
sponding to a we obtain the relation ()2.1|) . To prove the second part we express 
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Bi 
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A^KA 
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" K" 
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where the two partitioned matrices on the right hand side are the matrices a p'a 1 and 
up" a~ x . Now construct vectors t>j, w iy i = 1, 2, . . . , n such that 

ap'a- 1 = {((VilVj))), z,je{l,2,...,n} (2.4) 
ap"a- 1 = (((WilWj))), i,j 6 {1,2,. ..,n}. (2.5) 

Let |0 >, |1 > be the standard orthonormal basis of C 2 . Define 

\<Pi >= -p(|w< > |0 > +\wi > |1 >), 1 < i < n. (2.6) 
v2 

Then we have 

< ^il^i > = ^(( v i\ v j) + ( w i\ w j) 

= (ui\uj) for all i, 2, . . . , n}. 

Thus the correspondence itj — * ifi { is an isometry. Hence by (|2.3|) we have 

¥> k+j = a lj ip 1 + a 2j <p 2 + ■■■ + a kj cp k , 1 < j < n - k. 

Substituting for the (p^s from ()2.6|) and using the orthogonality of |0 > and |1 > we 
conclude that 

k 

\ v k+j > = y^jjlvj >, (2.7) 

i=l 
k 

\ w k+j > = ^2a ij \w i >. (2.8) 

i=l 

Putting 

= (((VilVj))), /../ c {1.2 /,•} 

AT" = (««;>,•))), i,je={l,2,...,fc} 

and substituting and (jHHJ) in (JO) and (Q we obtain B 1 = K'A, C x = A^K'A, 
B 2 = K"A, C 2 = A^K" A. Thus we have fl2HJ). ■ 
Let 7i\, 7i 2 be two complex Hilbert spaces of finite dimension d\,d 2 and equipped 
with orthonormal bases {e 1; e 2 , . . . , e dl }, {fi,f 2 ,---,fd 2 } respectively. Consider the 
tensor product 7i = Hi <S> H 2 equipped with the orthonormal basis = ej <8> fj with 
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the ordered pairs ij in the lexicographic order. For any operator X on H we associate 
its marginal operators Xi in Hi by putting 

Xi = Tr H2 X, X 2 = Tr ni X 

where Tr-^ stands for the relative trace over Hi. If p is a state on H, i.e., a positive 
operator of unit trace, then its marginal operators are states in Hi and H 2 . Now we fix 
two states pi and p 2 m 'Hi an d H 2 respectively and consider the compact convex set 

C(pi, P2) = {p\p a state 011H with marginals pi and p 2 in Hi and 7i 2 respectively. } 

in B(H). Let £(pi,p 2 ) C C(p±, p 2 ) be the set of all extreme points in C(pi,p 2 ). 

Proposition 2.2 Let p e £(pi,p 2 ). Then p is singular. 

Proof: Suppose p is nonsingular. Choose nonzero hermitian operators Li in Hi with 
zero trace. Then for all sufficiently small and positive e, the operators p ± eLi ® L 2 are 
positive definite. Since the marginal operators of L\ ® L 2 are 0, both of the operators 
p ± eLi ® L 2 belong to C(pi, p 2 ) and 

p = i ((p + eLx ® L 2 ) + (p - eLi ® L 2 )) 

and p is not extremal. ■ 



Proposition 2.3 Let n = did 2 , p E C(pi,p 2 ), rank p = k < n and let cr be a permu- 
tation of the ordered basis {fl^} of H such that 



1 _ 


K 


KA 






A^KA 



(2.9) 



where K is a strictly positive definite matrix of order k. Then, in order that p e £(pi, P2) 
it is necessary that there exists no nonzero hermitian matrix L of order A; such that both 
the marginal operators of 



1 


L 


LA 




A^L 


A^LA 



(2.10) 



vanish. 
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Proof: Suppose there exists a nonzero hermitian matrix L of order k such that both 
the marginals of the operator (j2.10j) vanish. Since K in (|2.9|) is nonsingular and positive 
definite it follows that for all sufficiently small and positive e, the matrices K ± e L are 
strictly positive definite. Hence 



1 


K + eL 


(K + eL) A 


a + a 1 


K-eL 


(K- 


eL)A 




A\K + eL) 


A\K + eL)A 


A\K-eL) 


A\K 


-eL)A 



a 



where each summand on the right hand side has the same marginal operators as p. 
Furthermore 



(K±eL) [I\A] > 0. 



K±eL 


(K±eL) 




' I 


A^(K±eL) 


A\K±eL)A 







Thus p is not extremal. 

Corollary Let p G p 2 )- Then rank p < \/df + d% — 1. 



Proof: Let rank p = k. By proposition 2.2, k < n. Since p is a positive definite matrix 
in the basis {g^} such that crpa^ 1 can be expressed in the form (|2.9|) . The extremality 
of p implies that there exists no nonzero hermitian matrix L of order k such that the 
matrix (j2.10j) has both its marginals equal to 0. The vanishing of both the marginals of 
(j2.1()|) is equivalent to 



Tr 



a 



1 


L 


LA 




. A]L 


A^LA 



a (Xi <g> / (2) + ® X 2 ) = 



(2.11) 



for all hermitian operators Xj in 7ij, JW being the identity operator in 7ij. Equation 
()2.11|) can be expressed as 



Tr L [I k \A\ a (X 1 <g> 1^ + /« <g> X 2 ) a' 1 



At 



0. 



In other words L is in the orthogonal complement of the real linear space 

h 



V = { [I k \A]a (Xi ® 1^ + J« ® X 2 ) a" 1 



A f 



Xj hermitian in 7ij, i = 1,2 
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with respect to the scalar product (L\M) = TrLM between any two hermitian matrices 
of order k. Thus the extremality of p implies that V 1 - = {0}. The real linear space of all 
hermitian matrices of order k has dimension k 2 . The real linear space of all hermitian 
operators of the form X l ® 1^ + ® X 2 is d\ + d 2 - 1 . Thus k 2 = dim£> < dj + d 2 ,-!. 



Proposition 2.4 Let p G C (pi, p 2 ) , k, a, K, A be as in Proposition 2.3. Suppose there 
is no nonzero hermitian matrix L of order k such that both the marginal operators of 



1 


L 


LA 






A^LA 



vanish. Then p G £(pi,p 2 ). 

Proof: Suppose p^8(pi, p 2 ). Then there exist two distinct states p',p" in C(pi,p 2 ) such 
that 

P=\(p' + p"), p'±p"- 

Since rank p — kit follows from Proposition 2.1 that there exist positive definite matrices 
K' , K" of order k such that 



ap^a 



1 _ 


K# 


K*A 




_ A^K* 


A^K*A 



where (p#, i*T#) stands for any of the three pairs (p, K), (p', K'), (p", K"). Since p' ^ p" 
and hence ap'a' 1 ^ ap"^ 1 it follows that K' ^ K" . Putting L = K' - K" ^ we 
obtain a nonzero hermitian matrix L of order k such that both the marginal operators 
of 



1 


L 


LA 






A^LA 



vanish. This is a contradicton. ■ 
Combining Proposition 2.3, its Corollary and Proposition 2.4 we have the following 
theorem. 

Theorem 2.5 Let Hi, H 2 be complex finite dimensional Hilbert spaces of dimension di, 
d 2 respectively. Suppose C(pi, p 2 ) is the convex set of all states p in H = Hi <E> H 2 whose 
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marginal states in Hi and H 2 are pi and p 2 respectively. Let {e^}, {fj} be orthonormal 
bases for Hi, H 2 respectively and let g tj = e^® f -, i — 1, 2, . . . , di, j = 1, 2, . . . , d 2 be 
the orthonormal basis of in the lexicographic ordering of the ordered pairs ij. In order 
that an element p in C(pi,p 2 ) be an extreme point it is necessary that its rank k does 
not exceed \/ d\ + d\ — 1. Let a be a permutation unitary operator in 7i, permuting the 
basis and satisfying 

apa~ 



1 _ 


if 


KA 






A^KA 



where K is a strictly positive definite matrix of order k. Then p is an extreme point of 
the convex set C(pi, p 2 ) if and only if the real linear space 



A t 



Xi hermit ian in Hi, i — 1, 2 



V = j [I fc |A]<7 (X! ® / (2) + ® X 2 ) ff- 1 

coincides with the space of all hermitian matrices of order k. 

Proof: Immediate from Proposition 2.3, its Corollary and Proposition 2.4 



3 The case Hi = H 2 = C 2 

We consider the orthonormal basis 



" 1 " 


, ll>= 


" " 











1 



in C 2 and write 

\xy >= \x > <S>\y > for all x, y G {0, 1}. 

Then e x = 1 00 >, e 2 = 1 01 >, e 3 = 1 10 >, e 4 = 1 1 1 > constitute an ordered orthonormal 
basis for C 2 <8> C 2 . For any state p in C 2 <S> C 2 define 

K p ((x,y),(x',y')) = (xy\p\x'y') x,y,x',y' G {0,1}. (3.1) 

If p has marginal states pi, p 2 then 

K p ((x,0),(x',0)) + K p ((x,l),(x',l)) = (x\ Pl \x'}, (3.2) 
K p ((0,y),(0,y')) + K p ((1, y), (1, y 1 )) = (y\p 2 \y') (3.3) 
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for all x, y, x', y' in {0, 1}. If p is an extreme point of the convex set C(pi, p^) it follows 
from Theorem 2.5 that the rank of p cannot exceed \fl. In other words, every extremal 
state p' in C(pi, p-i) has rank 1 or 2. When p\ = p 2 = \l we have the following theorem : 



Theorem 3.1 Let Ji\ = Ti.2 = C 2 . A state p in C(^I, is an extreme point if and only 
if p = \Q >< Q\ where 

|fl >= (|0 > ®|^o > +|1 > #i >) , 
{|-0o l^i >} being an ortho normal basis of C 2 . 

Proof: We shall first show that there is no extremal state p of rank 2 in C(|J, -/). 
To this end choose and fix a state p of rank 2 in C(\l, |J). Then the right hand sides 
of (|3.2|) and (|3.3|) coincide with hS xx > and respectively and in the ordered basis 

{ej, 1 < j < 4} the positive definite matrix if p of rank 2 in (j3.1j) assumes the form 



2 





y 


Z 


1-a 
2 


t 


-y 


t 


l-a 
2 


—x 


-y 


—X 


a 
2 



(3.4) 



for some < a < 1, x, y, z, t E C The fact X p has rank 2 implies that one of the following 
three cases holds : 



(1) 
(2) 



x 



x 

l-a 



y ¥ 



is strictly positive definite 



is strictly positive definite 



\y\ 2 = a ^ a ^ and one of the matrices 



2 Z 

™ a 
x 2 



(3) \x\ 2 

positive definite. 
We shall first show that case (3) is vacuous. We assume that 



l-a 
2 



l-a 
2 



is strictly 



I 1 2 I 1 2 

f = \y\ 



all — a) 



-, \z\ < — , rankK = 2. 
4 1 1 4 f 



(3.5) 



conjugation by the unitary permutation matrix corresponding to the permutation (1)(24)(3) 
brings (|3.4|) to the form 



a 
2 


z 


y 


X 


Z 


a 
2 


—x 


-y 


y 


—X 


l-a 
2 


i 


X 


-y 




l-a 
2 



(3.6) 



with rank 2. By Proposition 2.1 this implies that 

t 



l-q 
2 



t 



l-q. 
2 



A*KA 



where 



A = K- 1 


y x 


, K = 


a 
2 


z 




_-x -y _ 




Z 


a 

2 . 



2 c ' y 



^ e ty , substituting the expressions of ()3.8|1 in 



Putting x 

equating the 11-entry of the matrices on both sides of (3.7) we get 

2 



(3.7) 

(3.8) 
flSl|) and 



+ z e 



-i(0+<p) 







and therefore \z\ 2 = a contradiction. 

The case \t\ 2 < is dealt with in the same manner. 

Now we shall prove that p is not extremal. Express (|3.4j) as 



where 



K 



x 



X 
l-a 



K 


KA 


A^K 


A^KA 




= K- X 



A ] KA = dK~ 



d 



ail — a) 



y z 
t -y 



- \xr > o 



(3.9) 

(3.10) 
(3.11) 



This implies the existence of a unitary matrix U such that 



K?A = d*UK~*. 
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From ()3.10|) we have 



V Z =KA = d l l 2 K l l 2 UK- 1 ' 2 . 
t -y _ 

Hence Tr U = 0. Since U is a unitary matrix of zero trace it has the form 



U = e ld V 



where V is a selfadjoint unitary matrix of determinant —1. In particular 

A = d^e^K-^VK- 1 ' 2 

where V is selfadjoint and unitary. We now examine the linear space 

h 



V= { [I 2 \A] (X 1 ®I 2 + I 2 ®X 2 )[ 



A*- 



Xj is hermitian for each i 



(3.12) 



(3.13) 



In the ordered basis {ej, j = 1,2,3,4} it is easily verified that X 1 <g> I 2 + I 2 <E> X 2 in V 
varies over all matrices of the form 



X hermitian, p, q G R, r G 



" X + pI 2 


rl 2 


fl 2 


X + qI 2 



Thus 



V = { X + AXA ] + rA ] + fA + qAA ] +pl\X hermitian, p, qeR, r G C} 



We now search for a hermitian matrix L of order 2 in V x with respect to the scalar 
product (Xi\X 2 ) = Tr XiX 2 for any two hermitian matrices of order 2. In other words 
we search for a hermitian L satisfying 



Tr L = 0, Tr L K~ X I 2 V K 1 ' 2 = 

Tr L (X + dK-^VK-^XK-^VK' 1 / 2 ) = 



(3.14) 



for all hermitian X. (Here we have substituted for A from (|3.12jl ). 

Note that VdK^V K- 1 / 2 = B is a hermitian matrix of determinant —1. Thus 
(I3.14|) reduces to 

Tr L = 0, Ti LB = 0, L + BLB = 0. (3.15) 
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The matrix B can be expressed as 



B = WDW 1 

where W is unitary and 

r a 

-or 1 

Then for any £ 6 C the hermitian matrix 



a > 0. 



£ 



w 



satisfies (|3.15j) . In other words V 1 - ^ {0} and therefore the linear space T> in ()3.13j) is 
not the space of all hermitian matrices of order 2. Hence by Theorem 2.5, the state p is 
not extremal. 

Thus every extremal state p in |J) is of rank 1. Such an extremal state p has 

the form 

p = \Q >< n\ 



where 



|f2 > = ^2 a xy\xy >, 



1. 



The fact that \Q >< Q\ has its marginal operators equal to \l implies that ((a X y)) 
^(( u xy)) where {{u xy )) is a unitary matrix of order 2. Putting 



we see that 



y=0 



|fi>=- ? =(|0>|Vo>+|l>|^i >) 



(3.16) 



where {|0 >,|1 >} is the canonical orthonormal basis in C 2 and {l^o >> |"0i >} * s 
another orthonormal basis in C 2 (which may coincide with {|0 >, |1 >}). Varying the 
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orthonormal basis {\if)o >, l^i >} °f ^ 2 i n Q3-16|) we get all the extremal states of C(|/, \l) 
as \ fl >< 01. ■ 
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